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Abstract 

We derive the non-equilibrium single-particle momentum distribution function of a hadron reso- 
nance gas. We then study the effects that this newly derived expression can have in the freeze-out 
description of fluid-dynamical models of heavy ion collisions and compare it with the method 
traditionally employed, the 14-moment approximation. 



1. Introduction 

Fluid-dynamical models have been able to successfully describe the transverse momentum 
spectra and azimuthal transverse momentum anisotropics of particles observed in ultrarelativis- 
tic heavy ion collisions. One of the main ingredients of fluid-dynamical models is the so-called 
freeze-out procedure in which the fluid degrees of freedom are matched to kinetic or particle de- 
grees of freedom. Since experiments at RHIC and LHC measure particles and not fluid elements, 
the freeze-out is an essential step towards the comparison with experimental data. 

The Cooper-Frye formalism, usually applied to describe this matching process, demands the 
knowledge of microscopic information: the (single-particle) momentum distribution function of 
hadrons on the hypersurface in which freeze-out is performed (usually a constant temperature 
hypersurface). While this is well known for ideal fluids (Bose-Einstein or Fermi-Dirac distribu- 
tions in the local rest frame of the fluid), it has remained an open problem for viscous fluids. We 
remark that, so far, Israel and Stewart's (IS) simple ansatz for single component systems is still 
being employed in most fluid-dynamical calculations. 

In this case, the non-equilibrium correction to the momentum distribution of the /-th hadronic 
species , Sf^\ is assumed to be fT] 

<5/^'^ = ^uyl^kj, (1) 

where ^ is the four-momentum of the corresponding hadron and s^y is an expansion coefficient 
that should be matched to the fluid-dynamical variables. 

The above expression has two distinct approximations: 1) the momentum dependence was 
assumed to be quadratic and 2) the coefficient Sjjy was assumed to be the same for all hadronic 
species, being matched to the shear stress tensor. Here, e, P, and 

T are the energy density, thermodynamic pressure, and temperature of the fluid. In principle, both 
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these assumptions are incorrect. The limitations of the first assumption were already investigated 
in some works [2]. On the other hand, the second assumption and its domain of validity were 
just investigated in Ref. The purpose of this work is to further elaborate the studies regarding 
how (5/|['' actually depends on the different particle species. 

2. Method of moments 

We use the method of moments, as developed in Ref. fT], to compute the df^'' of a multi- 
component system without making any a priori assumption regarding its momentum dependence 
and the particle dependence of the expansion coefficients. First we factorize 6f^^ in the following 
way = /ok^k^k'' where f^^ is the local equilibrium distribution function, f^'^ = 1-1- af^^ 
(a = -1/1 for fermions/bosons), and (/ij^' is an out-of-equilibrium contribution. Next, we expand 

in terms of its moments using a complete and orthogonal basis constructed from particle four- 
momentum, A^, and fluid four-velocity, u'^. As done in Ref. f^], we use an expansion basis with 
two basic ingredients: 1) the irreducible tensors \,kY^ ,k'f kj\k\^ k'.kf\ - ■ ■ , analogous to the 
well-known set of spherical harmonics and constructed by the symmetrized traceless projection 
of A^' ■ ■ ■ A^'", i.e., kf' ■ ■ ■ kj'-^ = Kl'-''-v:k-' ■ ■ ■ k--, and 2) the orthonormal polynomials P^^^ = 
Zi"=o ^nr' ("/i^) ' which are equivalent to the associated Laguerre polynomials in the limit of 
massless, classical particles 13]. 

Then, the distribution function becomes, 

oo oo 

Jk ~ /ok ^ /0k/()k Zj Zj 'k ''f" ^ ' 

(=0 «=0 

where we introduced the energy-dependent coeflicients, T^^''"' = [A^f V^!] 2,"=,, a,™' ^'*k ^' 

The fields f^l^ ''' can be determined exactly using the orthogonaUty relations satisfied by the 

expansion basis and can be shown to correspond to the irreducible moments of 6f^\ 

f^:;" ^ {Ekkf" ■ ■ ■ kf)^ , = J dK, (. . .) (3) 

where gi is the degeneracy factor of the /-th hadron species and dKi - gicf^l [{2nf]^\ As 

long as this basis is complete, the above expansion fuUy describes , no matter how far from 
equilibrium the system is. 

Here, we are interested only on the effects arising from the shear-stress tensor. For this case, 
it is enough to fix = 2 (take only irreducible second-rank tensors) in the expansion above, 
i.e., neglect all scalar terms, e.g. bulk viscous pressure, irreducible first-rank tensors, e.g. heat 
flow, and tensors with rank higher than two (that never appear in fluid dynamics). The next 
approximation is the truncation of the expansion in momentum space, keeping only the term 
with n = (with t - 2 already fixed). Then, we obtain (for classical particles) 

fW - f® + J^^k- h (4) 

4 -4k + 2(£. + p.)r2^ '^'.^V- (V 

Above, TT^'' = p^Q, £,, and P,- are the shear-stress tensor, the energy density, and the thermody- 
namic pressure of the i-th particle species, respectively. Note that by keeping only the term with 
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n - (for { - 2)ws have the same momentum dependence obtained in the IS formula. However, 
our expansion coefficients are not independent of the particle type. In this formalism, this did 
not have to be assumed, but appeared naturally as a consequence of the orthogonality relations 
satisfied by the basis. 

In order to apply this expression to describe freeze-out, further approximations are required. 
This happens because in fluid dynamics we only evolve the total shear-stress tensor of the sys- 
tem (tt'"' - 2, 71^* ) and do not know, from fluid dynamics itself, how it divides into the individual 
shear-stress tensors of each hadron species (ti'!^). In order to extract this information, it is manda- 
tory to know how transient relativistic fluid dynamics emerges from the underlying microscopic 
theory. 

Basically, transient relativistic fluid dynamics is derived as the long, but not asymptotically 
long, time limit of the Boltzmann equation. In this limit, it is possible to relate the shear-stress 
tensor of individual particle species with the total shear stress tensor, n^"^, and the shear tensor, 
_ giM 11^)^ of the system. This was done in Ref. for single component systems, and can 
be easily extended to multi-component systems. Here, we just write down the solution and leave 
the detailed derivation of this relation to a future work. The solution is. 



lO 



lO 



cr^"- = aiTT"" + /3iCT^\ (5) 



where is the shear viscosity of the i-th species, while rj - Yiili is the total shear viscos- 
ity. The matrices Q are defined in such a way as to diagonalize the collision integral M, 
Dr^ - diag(^i, • ■ • ,Xn), with;i',- being the eigenvalues of M. With the truncation in momen- 
tum assumed above and considering only elastic 2-to-2 collisions between classical particles, M 
has the following simple form 

= J['6'< + C', (6) 
where y'^ = 1 - (1/2) 6'-' and W'^ , , , , is the transition rate. 



3. Simple model of hadrons 

In order to compute the coefficients a, and y8, appearing in Eq. |5]we must provide a set of 
hadronic cross sections. In this work, we estimate these coefficients for the first time using a 
simple hadronic model, in which all hadrons have the same constant cross section, cr,y - 30 mb. 
We consider only elastic collisions between the hadrons and include all hadrons up to a mass of 
1 .8 GeV. In this case, a, and /3, are actually independent of the value chosen for the cross section 
and all the difference between the various hadrons species originate solely from their different 
masses. In this simple scheme, we computed a, and y6, for all hadron species up to 1 .8 GeV and 
included all the decays of unstable resonances. 

In order to probe the effect this more sophisticated df^^ can have on heavy-ion observables, 
we computed the elliptic flow of pions, kaons, and protons using the fluid-dynamical model ap- 
phed in Ref. ISJ, with exactly the same parameters (GLmix initialization, HH-HQ shear viscosity 
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parametrization), but considering several choices of df^\ The resutls are shown in Figs.l. The 
different lines correspond to the following cases: 1) solid line uses Israel-Stewart's ansatz, 2) 
dotted line uses the 6/^'' derived in this paper, and 3) dashed line uses the 6f^^ derived in this 
paper in the Navier-Stokes limit (n^^ = Irjcr^'^). 
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Figure 1 : Differential elliptic flow for pions, kaons, and protons for three different centrality classes, (a) 20 - 30 %, (b) 
30 - 40 %, (c) 40 - 50 %, and various choices of c5/'''. 



We see that the particle dependence of the coefficients do not have a strong effect on the 
differential elliptic flow, which looks very similar to the one computed using Israel-Stewart's 
ansatz. At very non-central collisions, the difference is larger, but even then, cannot be considered 
as a crucial effect. Note that one of the reasons that makes this difference small is that we are 
not in the Navier-Stokes limit at freeze-out: the elliptic flow computed assuming the Navier- 
Stokes limit actually deviates more strongly from the one computed using Israel-Stewart's ansatz. 
Furthermore, note that there is a qualitative particle dependence in our result, the elliptic flow of 
pions is below the one computed with Israel-Stewart's ansatz while the elliptic flow of kaons and 
protons are always above. 

We remark that these results arise from an over simplified hadronic description and this in- 
sensitivity to the particle species might be coming just from the simple choice of hadronic cross 
sections assumed. In future calculations, we will come back to this discussion using more real- 
istic cross sections for the hadrons and also including inelastic colhsions between the hadrons as 
well. 

The work of H.N. was supported by the Academy of Finland, Project No. 133005. G.S.D 
thanks D. Molnar for helpful discussions. 
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